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ABSTRACT 

Context. Low-mass stars have been recognised as promising targets in the search for rocky, small planets with the potential of 
supporting life. As a consequence, Doppler search programmes using high-resolution spectrographs like HARPS or HARPS-N are 
providing huge quantities of optical spectra of M dwarfs. However, determining the stellar parameters of M dwarfs using optical 
spectra has proven to be challenging. 

Aims. We aim to calibrate empirical relationships to determine accurate stellar parameters for early M dwarfs (spectral types MO- 
M4.5) using the same spectra that are used for the radial velocity determinations, without the necessity of acquiring IR spectra or 
relying on atmospheric models and/or photometric calibrations. 

Methods. Our methodology consists in the use of ratios of pseudo equivalent widths of spectral features as a temperature diagnostic, 
a technique largely used in solar-type stars. Stars with effective temperatures obtained from interferometric estimates of their radii are 
used as calibrators. Empirical calibrations for the spectral type are also provided. Combinations of features and ratios of features are 
used to derive calibrations for the stellar metallicity. Our methods are then applied to a large sample of M dwarfs that are currently 
being observed in the framework of the HARPS GTO search for extrasolar planets. The derived temperatures and metallicities are 
used together with photometric estimates of mass, radius, and surface gravity to calibrate empirical relationships for these parameters. 
Results. A large list of spectral features in the optical spectra of early M dwarfs was identified. Erom this list the pseudo equivalent 
width of roughly 43% of the features shows a strong anticorrelation with the effective temperature. The correlation with the stellar 
metallicity is weaker. A total of 112 temperature sensitive ratios have been identified and calibrated over the range 3100-3950 K, 
providing effective temperatures with typical uncertainties of the order of 70 K. Eighty-two ratios of pseudo equivalent widths of 
features were calibrated to derive spectral types within 0.5 subtypes for stars with spectral types between K7V and M4.5V. Regarding 
stellar metallicity, 696 combinations of pseudo equivalent widths of individual features and temperature-sensitive ratios have been 
calibrated, over the metallicity range from -0.54 to -1-0.24 dex, with estimated uncertainties in the range of 0.07-0.10 dex. We provide 
our own empirical calibrations for stellar mass, radius, and surface gravity. These parameters are found to show a dependence on the 
stellar metallicity. Eor a given effective temperature, lower metallicities predict lower masses and radii, as well as larger gravities. 
Conclusions. 
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1. Introduction 

Ratios of equivalent widths or central depths of absorption lines 
with different excitation potentials have been widely used as 
temperature indicators in different kind of stars including solar- 
type (e.g. Gray & Johanson 1991; Gray 1994; Caccin et al. 2002; 
Kovtyukh et al. 2003; Biazzo et al. 2007; Montes et al. 2007; 
Sousa et al. 2010; Datson et al. 2012, 2014a,b), giant stars (e.g 
Gray 1989; Strassmeier & Schordan 2000), and supergiants (e.g. 
Kovtyukh 8c Gorlova 2000). To the best of our knowledge, this 
technique has however not been extended to the low-mass stars 
regime most likely due to the difficulties in analysing their opti¬ 
cal spectrum, mainly covered by molecular bands (in particular 
TiO and water) which blend or hide most of the atomic lines 
commonly used in the spectral analysis of solar-type stars. Fur¬ 
thermore, M dwarfs are intrinsically faint in the optical. 

The accurate determination of stellar parameters for M 
dwarfs has proven to be a difficult task. Regarding stellar metal- 
licity, some photometric calibrations based on optical and near- 
infrared photometry exist, a technique pioneered by Bonfils et al. 
(2005) and updated by Johnson & Apps (2009); Schlaufman 8c 
Laughlin (2010) and more recently by Neves et al. (2012, here¬ 
after NE12), although they require accurate parallaxes and mag¬ 
nitudes which are usually available only for nearby and bright 
stars. Another common technique to characterise M dwarf metal- 
licities is based on the use of spectroscopic indices which mea¬ 
sure the relative strength of the TiO molecular band with respect 
to the CaH molecular bands near 7000A (Lepine et al. 2007; 
Dhital et al. 2012; Lepine et al. 2013). Since the M dwarfs spec¬ 
tral energy distribution peaks at infrared wavelengths, some pre¬ 
vious works have performed a search for spectral features and 
indices in this spectral region. In particular, Rojas-Ayala et al. 
(2012) use the equivalent width of the Na i and Ca i triplet and 
the H 2 O-K 2 index in the K band of the spectra. This method¬ 
ology has been also applied by Terrien et al. (2012) and Mann 
et al. (2013a) providing calibrations for H and J/optical spectral 
bands, respectively. Large samples of M dwarfs have been char¬ 
acterised by means of near-infrared indices in the recent works 
by Newton et al. (2014a) and Gaidos et al. (2014). On the other 
hand, spectral synthesis has been tested in several works, usually 
on small number of stars focusing mainly on strong atomic lines 
or on spectral windows known to be less affected by molecular 
lines (Woolf & Wallerstein 2005; Bean et al. 2006; Maness et al. 
2007; Onehag et al. 2012; Rajpurohit et al. 2014). 

Concerning the effective temperature, very few M dwarfs are 
bright enough for a direct measurement of their radii (e.g. Boy- 
ajian et al. 2012, hereafter BO 12), a technique pioneered by Se- 
gransan et al. (2003). The most common technique for deter¬ 
mining the effective temperature of an M dwarf is the compar¬ 
ison of observed spectra with models atmosphere (e.g. Lepine 
et al. 2013; Gaidos et al. 2014). Casagrande et al. (2008, here¬ 
after CA08) provides optical/near-infrared photometric calibra¬ 
tions based on an extension of the infrared flux method (IRLM) 
for LGK dwarfs from Casagrande et al. (2006) to M dwarfs. 
However, significant systematic differences between tempera¬ 
tures based on CA08 calibrations and temperatures based on in¬ 
terferometric radii measurements have been recently noted by 
Mann et al. (2015). 


Telescopic Nazionale Galileo (TNG) Archive (programmes ID CAT- 
147, and A27CAT_83). 

** Our computational codes including the full and more 
detailed version of Tables 3, 4 and 6 are available at 
http://www.astropa.inaf.it/~jmaldonado/Msdlines.html 


Despite the intrinsic difficulties in their characterisation, low- 
mass stars are nowadays at the centre in the search for small, 
rocky planets with the potential capability of hosting life (e.g. 
Dressing 8c Charbonneau 2013; Sozzetti et al. 2013) In particu¬ 
lar, the radial velocity searches currently ongoing with HARPS 
at La Silla, and in the framework of the Global Architecture of 
Planetary Systems project* (GAPS; Covino et al. 2013) at the 
Telescopio Nazionale Galileo (TNG) with HARPS-N, are pro¬ 
ducing a large quantity of high resolution and high signal-to- 
noise ratio spectra. Lxoplanet searches would certainly benefit 
from a methodology to determine accurate stellar parameters us¬ 
ing the same spectra that are being used for the radial velocity de¬ 
terminations, i.e., without the necessity of observing at infrared 
facilities (usually from space) or relying on atmospheric models. 
Lollowing this line of reasoning, a methodology to characterise 
M dwarfs from high resolution optical spectra by using pseudo 
equivalent widths of features has been presented in a recent work 
by Neves et al. (2014, hereafter NL14). 

The idea of pseudo equivalent width can be further exploited 
in order to calibrate empirical relationships for M dwarfs. This 
is the goal of this paper, in which we present a large database of 
empirical calibrations of spectral features, ratios of features and 
combinations of both with the aim of deriving Tefj, spectral type, 
and metallicity, for early M dwarf stars by using optical HARPS 
and HARPS-N spectra (wavelength range 383-693 nm). Unlike 
NL14 we take as reference temperature scale the one provided by 
B012 and not CA08. We also apply our method to derive spec¬ 
tral types. Lurthermore, we provide relations for stellar masses, 
radii, and surface gravities so these quantities can be obtained 
without using parallaxes or photometry. We use our methods to 
characterise in an homogeneous and coherent way a sample of 
the M dwarfs which are currently being monitored in the HARPS 
GTO radial velocity programme (Bonfils et al. 2013). Late M 
stars are excluded from this study since exoplanet search around 
these stars is difficult at optical wavelengths. 

This paper is organised as follows. Section 2 describes the 
spectroscopic data used in this work. Section 3 describes our 
methodology and how empirical calibrations for the main stel¬ 
lar parameters are built. These calibrations are then applied to a 
large sample of stars and results are compared with other tech¬ 
niques in Section 4. The derived temperatures and metallicities 
are used to build empirical calibrations for stellar masses, radii, 
and gravities in Section 5. Our conclusions follow in Section 6. 


This work makes use of HARPS (Mayor et al. 2003) and 
HARPS-N (Cosentino et al. 2012) spectra mostly taken from 
archive. Specifically, the data is taken from: i) The LSO pipeline 
processed LEROS and HARPS archive ii) The ESO Science 
Data products Archive and Hi) The TNG Archive The cor¬ 
responding ESO and TNG programme IDs are listed in the foot¬ 
note to the paper title. In addition to the data from archive, some 
HARPS-N spectra have been provided by the GAPS team. 

The instrumental setup of HARPS and HARPS-N is almost 
identical so data from both spectrographs can be used together. 

The spectra cover the range 383-693 nm (HARPS-N), and 378- 
691 nm (HARPS). Both instruments provide a resolving power 

' http://www.oact.inaf.it/exoit/EXO-IT/Projects/Entries/2011/12/ 
27_GAPS.html 

^ http://archive.eso.org/wdb/wdb/eso/repro/form 

^ http://archive.eso.org/wdb/wdb/adp/phase3_spectral/form7phase3_collection 
=HARPS 

http://ia2.oats.inaf.it/archives/tng 


2. Spectroscopic data 
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of R~ 115000. The spectra are provided already reduced using 
ESO/HARPS-N standard calibration pipelines. Typical values of 
the signal-to-noise ratio are between 50 and 90 (measured at ~ 
5500 A). Wavelengths are on air. The spectra were corrected 
from radial velocity shifts by using the IRAF ^ task dopcor. For 
this purpose, we used the accurate radial velocities (measured 
by applying the cross-correlation technique) which are provided 
with the reduced spectra. 

3. Methodology 

The optical spectrum of an M dwarf is a forest of lines and 
molecular bands usually heavily blended in which identifying 
individual lines (of moderate strength) or measure equivalent 
widths is a difficult task. In order to overcome this limitation 
we follow the idea depicted by NE14 and instead of considering 
equivalent widths or depths of lines, we consider pseudo equiv¬ 
alent widths (hereafter EWs) of features. A feature can be a line 
or a blend of lines. The pseudo equivalent width is defined as 
the traditional equivalent width, with the difference that it is not 
measured with respect to a continuum normalised to unit, but to 
the value of the flux between the peaks of the feature at each 
wavelength; 

EW (1) 

r pp 

where F.i is the stellar flux, Fpp denotes the value of flux between 
the peaks of the feature at each integration step, and AT is the 
spectra wavelength’s step. An estimate of the uncertainty in the 
measured EWs is given by: 


- — - -O-F,, ( 2 ) 

^ E PP ^ 

where < Epp > is the mean value of the flux between the peaks of 
the feature, and its corresponding standard deviation. Fig¬ 
ure 1 illustrates how the EWs are measured. 

An initial list of 4224 features was built taking as reference 
the spectra of the star G1 49, a M1.5V low-mass dwarf. Spec¬ 
troscopic observations of this star were carried out within the 
framework of the GAPS project with HARPS-N. Only the red 
region of the spectra (5300 - 6900A) was considered since the 
blue part of the optical spectrum of an M dwarf usually suf¬ 
fers from lower signal-to-noise ratio. Regions of the spectra 
affected by chromospheric activity and atmospheric absorption 
were avoided. 

Figure 2 shows the effects of effective temperature and 
metallicity on the EWs measurement. In the left panel a portion 
of the stellar spectra is shown for stars with similar metallici- 
ties but different Tefj. The same portion of the spectra is shown 
in the right panel, this time for stars of similar Tefr but different 
metallicities. 

A list of calibrators was built for each of the basic stellar 
parameters considered in the present work (Teff, spectral type, 
and metallicity), following different criteria as explained in the 
next subsections. 


^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Foundation. 



Wavelength (A) 

Fig. 1. Example of the measurement of EWs. The figure shows a portion 
of the spectrum of the star G1 49. The flux between the peaks of the 
features is measured over the blue lines. 


3.1. Effective temperature 

The accuracy of the Teff derived from temperature sensitive ra¬ 
tios is intimately tied to the accuracy of the temperature of the 
stars used as calibrators. We used as calibrators the sample of 
early M dwarfs with angular sizes obtained with long-baseline 
interferometry to better than 5% given in B012. 

BO 12 list 22 low-mass stars with spectral-types equal to or 
cooler than K5V, spanning a range of Teff between 3000 and 
4500 K. HARPS spectra were obtained for seven stars from the 
ESO archive, whilst HARPS-N spectra were obtained for three 
stars from the TNG archive. HARPS-N spectra for two further 
stars have been provided by the GAPS team. To these stars, we 
added two more from the recent work by von Braun et al. (2014) 
who analyze their stars in the same way as BO 12 

Mann et al. (2013b) revised the temperature scale of BO 12 
by noticing an issue regarding the determination of the bolomet- 
ric flux of the stars. Although the temperature differences are 
relatively small, we use the set of updated temperatures. For the 
two stars taken from von Braun et al. (2014), updated Teff values 
computed in the same way as in Mann et al. (2013b) are pro¬ 
vided in Newton et al. (2014b). The final list of Teff calibrators 
amounts to 14 stars whose parameters are listed in Table 1. 

Starting from our initial list of 4224 identified features, the 
EW of all features were measured in all calibration stars. In order 
to avoid possible dependencies on microturbulence, rotation, or 
stellar metallicity, we rejected features with EW < 20 mA or 
EW >120 mA in any of the calibration stars, thus excluding 
too weak and too strong features. We also rejected lines with 
relative errors i^-) larger than 2%. For every possible ratio of 
features a Spearman’s correlation test was computed to check 
whether the ratio is temperature sensitive or not. The ratios were 
selected with the only condition that the features are separated by 
no more than 15 A. This limit was set in order to avoid problems 
with scattered light correction or continuum normalisation. All 


® The authors list three M dwarfs but we were unable to obtain a 
HARPS or HARPS-N spectra of G1 649. 
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Wavelength (A) Wavelength (A) 

Fig. 2. Left: Portion of the stellar spectrum of several Tcff calibrators (see text for details) with similar metallicity (from -0.07 to -l- 0.15 dex) but 
different effective temperature: 3850 K (purple), 3701 K (blue), 3646 K (green), 3520 K (orange), and 3176 K (red). Right: Portion of the stellar 
spectrum of several metallicity calibrators (see text for details) with similar effective temperature (3450 - 3550 K) but different metallicity: -t-0.06 
dex (purple), -l-O.Ol dex (blue), -0.11 dex (green), -0.20 dex (orange), and -0.37 dex (red). For the sake of clarity, an offset of -0.40 in flux was 
applied between the spectra. 



_^_I_ I _^_ I _I_^_I_ 

3200 3400 3600 3800 4000 



3200 3400 3600 3800 4000 


Tetf (K) (K) 

Fig. 3. Examples of ratios of some features identified to be sensitive to Teff in early-M dwarfs. Stars are plotted using different colours according 
to their metallicity (using 0.10 dex length bins). Median uncertainties are shown in the left upper corner of each plot. The corresponding fits are 
also plotted. The features’ central wavelengths as well as the rms standard deviation of the residuals are given in each plot. 


ratios with a probability of correlation by chance lower than 2% 
were considered for further study ^. 

’ Ideally, for a ratio of lines to be temperature-sensitive the excitation 
potential of the lines, must differ as much as possible. This is because 
the EWs of lines with higher change with Tefj faster than those of lines 
with lowers values (Gray 1994). 


Following Kovtyukh et al. (2003) for each considered EW 
ratio we fitted the Teff-ratio relationship to several functions: a 
Hoerl function (Tefj = ab‘xr‘^), a modified Hoerl function (Tefr 
= ab'^’^xr'^), a power-law function (Teff = axr*’), an exponential 
law function Teff = axb‘, and a logarithmic function Teff = a H- 
bxln(r), where r = EW 1 /EW 2 is the ratio between the EW of two 
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Table 1. The effective temperature calibration sample. 


Table 3. Spectral Type calibration sample. 


Star 

SpT 

Teff 

(K) 

GJ338A 

MOV 

3953 +41 

GJ205 

M1.5V 

3850 ± 22 

GJ880 

M1.5V 

3731 + 16 

GJ176 

M2.5V 

3701 + 90 

GJ887 

M0.5V 

3695 + 35 

GJ526 

M1.5V 

3646 ± 34 

GJ15A 

M1.5V 

3602 + 13 

GJ412A 

MIV 

3537 + 41 

GJ436 

M3V 

3520 + 66 

GJ581 

M2.5V 

3487 + 62 

GJ725A 

M3V 

3417+ 17 

GJ699 

M4V 

3238 + 11 

GJ876 

M5V 

3176 + 20 

GJ725B 

M3.5V 

3142 + 29 


Notes. Effective temperatures are from Mann et al. (2013b), and spec¬ 
tral types for B012. For GJ176 and GJ876, Tefi values are from New¬ 
ton et al. (2014b), and spectral types from the GJ catalogue (Gliese & 
JahreiB 1991). 


GJ 

SpT 

Reft 

GJ 

SpT 

Reft 

185 

K7V 

h 

408 

M2V 

h 

686 

MOV 

h 

250B 

M2.5V 

k 

701 

MOV 

h 

581 

M2.5V 

h 

846 

M0.5V 

k 

352A 

M3V 

k 

720A 

M0.5V 

k 

436 

M3V 

k 

229 

MIV 

k 

752A 

M3V 

k 

412A 

MIV 

h 

569A 

M3V 

k 

514 

MIV 

h 

273 

M3.5V 

k 

570B 

MIV 

h 

643 

M3.5V 

k 

908 

MIV 

h 

734B 

M3.5V 

k 

205 

M1.5V 

k 

213 

M4V 

k 

15A 

M1.5V 

h 

402 

M4V 

k 

526 

M1.5V 

h 

699 

M4V 

k 

625 

M1.5V 

h 

83.1 

M4.5 

k 

220 

M2V 

k 

166C 

M4.5 

k 

382 

M2V 

k 

234A 

M4.5 

k 

393 

M2V 

h 





Notes. ^ k: standard from Kirkpatrick et al. (1991); h: standard from 
Henry et al. (1994). 


Table 2. Coefficients of our feature ratio-temperature relations. 
Columns (1) to (4) provide information about the features involved in 
the ratio (central wavelength and width in A), while columns (5) to (9) 
show the coefficients of the best-fitting relationships, their functional 
form, and the corresponding standard deviation of the Teff calibration. 
Only four examples are shown here. The same examples are shown in 
Figure 3. 


ii 

(1) 

A2i 

(2) 

/I? 

(3) 

A/i? 

(4)' 

a 

(5) 

b 

(6) 

c 

(7) 

func.^ 

(8) 

o-(K) 

(9) 

6435.00 

0.38 

6432.01 

0.34 

4529.63 

0.73 

-0.015 

MH 

49 

5720.46 

0.33 

5718.32 

0.28 

4842.67 

0.66 

0.89 

H 

53 

6226.67 

0.35 

6213.44 

0.29 

4350.09 

0.86 

- 

E 

54 

5708.61 

0.30 

5718.32 

0.28 

5122.92 

0.634 

0.92 

H 

54 


H: Hoerl; MH: modified Hoerl; PL: Power-law; E: Exponential; L: Logarithmic 


features. All fits were performed using a nonlinear least-squares 
fitting routine in IDL (MPFIT; Markwardt 2009) taking into ac¬ 
count the uncertainties in Teff. For each calibration we selected 
the function that produced the smallest standard deviation, re¬ 
taining only those calibrations with a standard deviation smaller 
than 75 K. The number of selected temperature-sensitive ratios 
amounts to 112. 

Given our relatively low number of calibrators we performed 
an additional check to ensure that the selected ratios are not cor¬ 
related with Teff simply by coincidence. We created 1000 series 
of simulated random effective temperatures and errors, keeping 
the media and the standard deviation of the original data. For 
each series of simulated data we repeted our analysis and com¬ 
puted the number of “suitable” calibrations. The results show 
that in 98% of the simulations our methodology does not recover 
any suitable Teff-calibration, whilst only in 0.8% of the simula¬ 
tions the number of obtained calibrations is larger than 10. We 
therefore conclude that it is very unlikely that our obtained Tgff- 
ratios are correlated with Teff just by chance. 

Some examples of the selected temperature-sensitive ratios 
are shown in Figure 3, whilst full details regarding the calibra¬ 
tions for the same examples can be found in Table 2. 


Table 4. Coefficients of our spectral type calibrations. Columns (1) and 
(2) show the wavelength of the corresponding features, columns (3) to 
(6) the coefficients of the fit, while column (7) gives the standard devia¬ 
tion of the calibration of spectral types. Only four examples are shown 
here. 


M 

(1) 

^2 

(2) 

ao 

(3) 

ai 

(4) 

a2 

(5) 

a3 

(6) 

(T 

(7) 

5468.06 

5481.92 

-3.88 

10.79 

-4.25 

0.48 

0.49 

5627.63 

5667.53 

5.11 

-0.75 

-2.18 

0.50 

0.41 

5467.06 

5481.92 

-3.80 

8.18 

-0.56 

-0.77 

0.48 

5467.06 

5512.54 

-5.09 

10.67 

-4.13 

0.53 

0.45 


3.2. Spectral Types 

A similar approach was followed to derive spectral types. 
HARPS or HARPS-N spectra were obtained for a sample of 33 
stars with homogeneously derived spectral types in Kirkpatrick 
et al. (1991) and Henry et al. (1994). The sample contains stars 
with spectral types between K7V and M4.5V. These stars are 
listed in Table 3. 

Spectral type-sensitive ratios were identified by means of a 
Spearman’s correlation test. For each of them a third order poly¬ 
nomial fit was performed between the numerical spectral-type 
index (with value 0.0 for MO; 0.5 for MO.5 and so on) and the 
ratio. A negative value implies that the star is a late-K dwarf in¬ 
stead of an M star being the index value for K7 equal to -1.0 (K7 
is the subtype preceding MO). 

The third order polynomial fit was preferred amongst other 
functional fits since we found it to give the lowest rms standard 
deviation. Our “final” selection of spectral type-sensitive ratios 
includes 82 ratios of features with a standard deviation lower 
than 0.5 spectral subtypes. The derived mean numerical spectral 
types are rounded to the nearest half integer. A couple of exam¬ 
ples are shown in Figure 4. Full details for some examples are 
given in Table 4. As for T^ff, we performed 1000 simulations 
with random spectral-type values. The results show that we are 
not able to find any suitable calibration in any of the simulations 
when random spectral-types are used. 
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Fig. 4. Spectral type as a function of two different spectral-type sensi¬ 
tive ratios. Stars are plotted using different colours according to their 
metallicity (using 0.10 dex length bins). A third order polynomial fit is 
shown. The features’ central wavelengths as well as the rms standard 
deviation on spectral type of the residuals are given in each plot. 


3.3. Metallicity 

A common approach to find metallicity calibrators for low-mass 
stars relies on the search of M dwarfs in common proper mo¬ 
tion pairs orbiting around a solar-type star with accurate spectro¬ 
scopic metallicity determinations. The basic assumption is that 
both stars are coeval and born in the same protostellar cloud so 
the metallicity of the secondary M dwarf is the same as the one 
of the primary star (e.g. Bonfils et al. 2005). Until very recently, 
only the most nearby and bright M stars have been searched for 
planets by means of the Doppler technique. As a consequence 
there is a lack of HARPS and HARPS-N spectra for most of the 
identified M dwarfs in binary systems around solar-type stars. 
To overcome this difficulty we built a list of 47 metallicity cal¬ 
ibrators with available HARPS spectra, known parallaxes, and 
magnitudes by using the most recent photometric calibration 
provided in NE12. This calibration is in turn based on metal¬ 
licity determinations from FGK primaries with an M dwarf sec¬ 
ondary. The sample of metallicity calibrators covers a wide range 
in metallicity from -0.54 to h- 0.24 dex with typical error bars of 
the order of < 0.05 dex. These stars are listed in Table 5. We 
caution that the uncertainties reported in Table 5 do not take into 
account the scatter in the NE12 calibration, which is of the order 
of-0.17 dex. 

The effects of metallicity on the EW of the features are en¬ 
tangled with the effects of Tefj, with Tefj as the primary driver 
of changes in the EW. This can be easily seen in the histograms 
in Figure 5. They show the distribution of the Spearman’s rank 
correlation factor of the EW with the stellar metallicity and with 
Teff. The figure shows that a significant fraction of the features, 
- 43%, shows a high (Spearman’s correlation factor < -0.80) 
anticorrelation with Tes, while only a relatively small fraction 
(~ 3%) shows a significant positive correlation. The correlation 
between stellar metallicity and EWs is generally less significant. 


Table 5. Metallicity calibration sample. [Fe/H] are computed from V 
magnitudes and parallaxes taken from the compilation of Bonfils et al. 
(2013) and 2MASS magnitudes (Cutri et al. 2003). The photometric 
calibration by NE12 is used. Errors are computed by propagating the 
uncertainties in the parallaxes and the photometry. 


Star 

[Fe/H] 

Star 

[Fe/H] 

GJl 

-0.37 + 0.04 

GJ551 

0.13 + 0.04 

GJ105B 

-0.13 + 0.03 

GJ555 

0.10 + 0.04 

GJl 76 

0.02 ± 0.04 

GJ569A 

0.02 + 0.03 

GJ205 

-0.03 + 0.19 

GJ588 

0.05 + 0.03 

GJ2066 

-0.10 + 0.03 

GJ618A 

-0.06 + 0.04 

GJ213 

-0.24 + 0.04 

GJ628 

-0.05 + 0.03 

GJ229 

-0.02 + 0.17 

GJ674 

-0.20 + 0.03 

GJ250B 

-0.10 + 0.04 

GJ678.1A 

-0.14 + 0.04 

GJ273 

-0.11+0.03 

GJ680 

-0.05 + 0.04 

GJ300 

0.06 + 0.03 

GJ682 

0.10 + 0.03 

GJ357 

-0.32 + 0.03 

GJ686 

-0.30 + 0.03 

GJ358 

0.05 + 0.03 

GJ693 

-0.29 + 0.03 

GJ367 

-0.05 + 0.04 

GJ701 

-0.19 + 0.03 

GJ382 

0.05 + 0.03 

GJ752A 

0.01 + 0.03 

GJ393 

-0.11 +0.04 

GJ832 

-0.17 + 0.04 

GJ413.1 

-0.11 +0.04 

GJ846 

-0.07 + 0.04 

GJ438 

-0.51+0.07 

GJ849 

0.24 + 0.04 

GJ447 

-0.26 + 0.03 

GJ876 

0.13 + 0.03 

GJ465 

-0.54 + 0.04 

GJ877 

-0.02 + 0.03 

GJ479 

0.05 + 0.04 

GJ887 

-0.35 + 0.14 

GJ514 

-0.11 +0.04 

GJ908 

-0.39 + 0.03 

GJ526 

-0.16 + 0.03 

HIP31293 

0.03 + 0.04 

GJ536 

-0.15 + 0.04 

LTT9759 

0.14 + 0.04 

GJ54.1 

-0.47 + 0.05 




The distribution of Spearman’s coefficients for metallicity shows 
a clear peak at -1-0.25 which drops almost to zero at -1-0.50, while 
at negative values it has a smooth tail down to -0.80. Effects of 
metallicity and effective temperature should therefore be consid¬ 
ered simultaneously. 

We searched for empirical relationships for metallicity as a 
function of features and ratios of features (i.e., an indicator of 
temperature) with the following analytical form; 

[Fe/H] = (A X EW) + iBxr) + C (3) 

where EW is the EW of a feature showing a strong-metallicity 
correlation, r is a temperature-sensitive ratio of features, and A, 
B, C are independent coefficients. We considered every combi¬ 
nation of features and ratios satisfying the condition that the cor¬ 
relation of EW with metallicity, and the correlation of r with Teff 
show at least a 98% of significance. Our final selection consists 
of 696 calibrations with standard deviation values between 0.07 
and -0.10 dex. We point out that these uncertainties should be 
considered as lower limits since they do not take into account 
possible systematic errors in the underlying NE12 calibration. 

As before, we performed a series of simulations using ran¬ 
dom metallicities and errors. In 84% of the simulations we do not 
find any suitable metallicity calibration, although in 7.5% of the 
cases the simulation finds a large number of calibrations (larger 
than 348). Some examples of our obtained metallicity calibra¬ 
tions are provided in Table 6. 

4. Comparison with other methods 

Our calibrations have been applied to a sample of 53 M dwarfs 
from the HARPS GTO M dwarf sample (Bonfils et al. 2013) for 
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Spearmon correlation foctor 

Fig. 5. Spearman’s correlation factor distribution of the EW with the 
stellar metallicity (blue histogram) and with the Tefr (red histogram). 
The sample of metallicity calibrators and the initial list of 4224 features 
is considered. Metallicities are computed using NE12 whilst effective 
temperatures are obtained with our methodology. 

Table 6. Coefficients of our metallicity calibrations. Column (1) 
shows the wavelength of the corresponding feature, column (2) the 
temperature-sensitive ratio, columns (3) to (5) the coefficients A, B, and 
C, while column (6) gives the standard deviation of the calibration. Only 
five examples are shown here. 


-tl 

(1) 

F/b 

(2) 

A 

(3) 

B 

(4) 

C 

(5) 

(T (dex) 

(6) 

6785.77 

6785.38/6799.25 

-0.0354 

-1.876 

1.753 

0.07 

6785.77 

6799.25/6785.38 

-0.0349 

0.227 

0.416 

0.07 

6785.77 

6785.38/6790.93 

-0.0359 

-0.930 

1.633 

0.07 

6785.77 

6785.38/6788.76 

-0.0375 

-1.163 

1.581 

0.07 

6785.77 

6790.93/6785.38 

-0.0356 

0.263 

0.613 

0.07 


which HARPS data have been obtained from the ESO Science 
Data Products Archive^. Only spectra with a median signal-to- 
noise ratio of at least 25 were considered. For stars with more 
than one spectrum available we took the one with the highest 
signal-to-noise ratio. No further restrictions were applied. The 
sample is composed of nearby (distance <11 pc), bright (V < 
12, Ks < 7), early-type M dwarfs (spectral types M0V-M4.5V). 
Our methods were applied to compute effective temperatures, 
stellar metallicities and to derive spectral-types. “Final” values 
for these parameters are the mean of the individual values from 
all the calibrations. All these quantities are provided in Table 8, 
which is available in the online version of this paper. 

Our results are compared with: i) A photometric scale, 
namely CA08 for Tgff and NE12 for [Fe/H]; ii) The recent work 
by Gaidos et al. (2014, hereafter GA14); and Hi) The values ob¬ 
tained with the methodology developed by NE14. 

4.1. Comparison of effective temperatures 

Photometric effective temperatures are derived from V magni¬ 
tudes from the compilation of Bonfils et al. (2013) and 2MASS 
(Cutri et al. 2003) photometry using the calibration provided by 
CA08. Computed errors take into account the propagation of the 
uncertainties of the 2MASS magnitudes as well as the accuracy 

* archive.eso.org/wdb/wdb/adp/phase3_spectral/form?phase3_collection 
=HARPS 


of the CA08 calibrations. The comparison between the photo¬ 
metrically derived temperatures and our spectroscopic ones is 
illustrated in Figure 6. There is a clear offset between our spec¬ 
troscopic estimates and the photometric temperatures, being the 
latter cooler than ours (the median difference ATeff = - T^^™ 

is -198 K with a rms standard deviation of 176 K). Our temper¬ 
atures can be converted into the CA08 scale by a linear trans¬ 
formation: Teff[CA08 scale] = (1.29 + 0.02)xTeff[this work] h- 
(-1271 + 89) K (dashed grey line). 

The reason for this discrepancy is the choice of Mann et al. 
(2013b) temperatures as our reference temperature scale. In or¬ 
der to test this, photometric Tefr values were computed for our 
sample of temperature calibrators (Table 1) following CA08. The 
comparison between Mann et al. (2013b) and CA08 tempera¬ 
tures is shown in Figure 7. It can be seen from this figure that 
CA08 temperatures tend to be systematically lower than those 
provided by Mann et al. (2013b). 

The discrepancy between CA08 values and interferometry- 
based temperatures has also been noted in a recent work by 
Mann et al. (2015). The difference ATgff between CA08 and 
interferometric-based temperatures noted by these authors is 160 
K (the CA08 temperatures being cooler) in agreement with our 
results. CA08 temperatures are obtained by extending the IRFM 
for FGK dwarfs from Casagrande et al. (2006) to M dwarfs. One 
of the assumptions of the IRFM is that a star can be approxi¬ 
mated as a blackbody for wavelengths beyond » 2 pm. Mann 
et al. (2015) argue that whilst this assumption is reasonable for 
warmer stars, it does not suit M dwarfs, which have significantly 
more flux in the near-infrared than predicted by a blackbody. 
As a consequence CA08 temperatures tend to be systematically 
lower, with increasing disparity at cooler temperatures where 
stars deviate more from the blackbody emission (see Figures 6 
and 7). Mann et al. (2015) also note that the temperature scale 
of the old version of the PHOENIX models used in CA08 differs 
from interferometric-based temperatures. 

We have also compared our temperatures with the values 
given by NE14. Since NE14 method is calibrated using the 
CA08 photometric relationship, the comparison of our temper¬ 
atures with NE14 shows results similar to the comparison with 
CA08 (Figure 6). 

Our sample contains 51 stars in common with the sample 
of GAM who determine effective temperatures by comparing 
low-medium resolution spectra with the BT-SETTF version of 
the PHOENIX model atmospheres (Allard et al. 2012a,b). Their 
procedure was calibrated using the stars listed in BO 12, although 
with the stellar bolometric fluxes computed as in Mann et al. 
(2013b). As can be seen in Figure 6, GAM temperatures tend to 
be slightly hotter than ours, specially for Tefj > 3400 K. For the 
coolest dwarfs in this sample, GAM temperatures depart from 
ours and tend to be smaller than ours. 

4.2. Comparison of metaiiicities 

We also compare our metallicities with those reported previ¬ 
ously in the literature. Values for the comparison are taken from 
the photometric calibration by NE12; from GAM who deter¬ 
mine metallicities following the method of Mann et al. (2013b) 
based on empirical calibrations between the strength of atomic 
and molecular spectroscopic features and stellar metallicity; and 
from NEM. The comparison is shown in Figure 8. 

The comparison reveals an overall good agreement between 
our metallicity estimates and those by NEM, GAM, and NEM. 
The median differences with these works are consistent with 
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Fig. 6. Teff values from the literature estimates versus the values ob¬ 
tained in this work. The upper panel shows the differences between the 
temperatures given in the literature and the values derived in this work. 
Median uncertainties in the derived temperatures are also shown. The 
symbol <> in the legend represents the median difference. The black 
continuous line represents the 1:1 relation whilst the grey dashed one 
represents the best linear fit between our estimates and those obtained 
using CA08 relationship (see text in Section 4.1). 



0) 


- 0.5 



[Fe/H] (This work) 


Fig. 8. [Fe/H] values from the literature estimates versus the values 
obtained in this work. The upper panel shows the differences between 
the metallicities given in the literature and those derived in this work. 
Median uncertainties in the derived metallicities are also shown. The 
symbol <> in tbe legend represents the median difference. The black 
continuous line represents the 1:1 relation whilst the grey dashed one 
represents the best linear fit between our estimates and those obtained 
using NE12 relationship (see text in Section 4.2). 



(K) Mann 

Fig. 7. Teff estimates based on the calibration provided by CA08, versus 
values from Mann et al. (2013b) and Newton et al. (2014b). The upper 
panel shows the differences. The symbol <> in the legend represents the 
median difference. The black continuous line represents the 1; 1 relation 
whilst the grey dashed one represents the best linear fit. 


zero and the scatter although somewhat large is consistent within 
the (also large) error bars. A linear fit between our metallici¬ 
ties and those obtained using NE12 relationship provides a slope 
slightly larger than one and a small difference in the zero point; 
[Fe/H][NE12 scale] = (1.22 + 0.04)x[Ee/H][this work] H- (0.02 
+ 0.01) dex (dashed grey line). 


4.3. Comparison of spectral types 

We finally compare the spectral types derived by us with those 
obtained by using the automatic procedure of the HAMMER 
spectral code (Covey et al. 2007, hereafter CO07). The CO07 
code was designed to classify stars in the Sloan Digital Sky Sur¬ 
vey Spectroscopic database, therefore before using it our spec¬ 
tra were degraded to a resolution R ~ 2000 by convolving them 
with a gaussian profile. We also caution that roughly half of the 
spectral-type sensitive band indices defined in CO07 are outside 
the HARPS spectral coverage. The comparison is shown in Fig¬ 
ure 9. It can be seen that there seems to be no significant differ¬ 
ences between our estimates and those by CO07, with an overall 
good agreement within + 1 spectral subtype (dashed lines in Fig¬ 
ure 9). West et al. (2011) and Lepine et al. (2013) found the auto¬ 
matic spectral types given by the Hammer code to be about one 
subtype earlier than the manual classification. While this effect 
is not evident in our comparison it can not be ruled out either. 

Figure 10 shows our derived effective temperatures as a func¬ 
tion of the spectral type. For comparison data from Kenyon & 
Hartmann (1995, Table A.5) is overplotted (red circles). The me¬ 
dian Teff-spectral type sequence from Fepine et al. (2013) is also 
shown (green squares). It can be seen that except for the pres¬ 
ence of some outliers, effective temperatures and spectral types 
are well correlated following the expected trend. In other words, 
our spectral types appear to be consistent with our temperature 
scale. Unlike Fepine et al. (2013), our data do not suggest the 
presence of a Tgff plateau in the spectral range Ml-M3 although 
we note that our sample is relatively small in comparison with 
the one in Fepine et al. (2013). 

We conclude that our metallicities and spectral types agree 
reasonably well with other literature estimates. Regarding ef¬ 
fective temperatures, there is a clear offset between BO 12 and 
CA08 scales, as explained. In summary, our methodology can 
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SpType (This work) 

Fig. 9. M subtype values obtained by using the HAMMER code versus 
the values obtained in this work. The upper panel shows the differences 
with the values given in the literature. Random values between ±0.2 
have been added to the Hammer values to help in the comparison. The 
symbol <> in the legend represents the median difference. The black 
confinuous line represents the 1:1 relation whilst the dashed ones corre¬ 
spond to ± 1 spectral subtype. 



Spectral Type 

Fig. 10. Effective temperature as a function of the spectral type. Eor 
clarity the spectral types are not rounded. A second order polynomial fit 
is shown. Possible outliers are removed by using a 2.5cr clipping pro¬ 
cedure. The coefficients of the fit as well as the rms standard deviation 
are given in the plot. Data from Kenyon & Hartmann (1995, Table A.5) 
is overplotted using red filled circles, whilst green squares represent the 
median Tefj-spectral type sequence from Lepine et al. (2013). 


be confidently used to characterise large samples of stars in an 
homogeneous way. 


5. Empirical relationships for stellar mass, radius, 
and gravity 

We made use of the temperature and metallicity values derived 
with our method to search for empirical relationships with the 
stellar evolutionary parameters namely, stellar mass, radius, and 
surface gravity. 

We derived our own mass-radius relationship by combining 
the stars with known interferometric radius from BO 12 and von 
Braun et al. (2014) with data from low-mass eclipsing binaries 
provided in the compilation by Hartman et al. (2014, Table 5). 
A 3cr clipping procedure was used to remove potential outliers. 
Our derived calibration is as follows: 


R = 0.0753 -H 0.7009 x M + 0.2356 x (4) 

where radius and masses are given in solar units and the rms 
standard deviation of the calibration is 0.02 Rq. The radius-mass 
plane is shown in Figure 11. 



Fig. 11. Stellar radius as a function of the stellar mass. The sample in¬ 
cludes stars with interferometric measurements of their radii from B012 
and von Braun et al. (2014) as well as low-mass eclipsing binaries from 
Hartman et al. (2014, Table 5). The best fit is also shown. A 3cr clipping 
procedure was used to remove outliers. The upper panel shows the dif¬ 
ferences between the radius derived with our fit and the radius given in 
the literature. Median errors in radius (not shown in the plot) are of the 
order of 0.006 R©. 

Values of stellar masses were obtained for each of our target 
stars following the relations based on near infrared photometry 
by Henry & McCarthy (1993). We chose this calibration since 
it is the same used by B012. These calibrations are provided 
in the CIT photometric system therefore, before applying them, 
2MASS magnitudes were converted into CIT magnitudes fol¬ 
lowing the transformations provided by Carpenter (2001). Once 
the stellar masses were computed, values of the radius were de¬ 
rived using Equation 4. Finally, surface gravities, log g, were de¬ 
rived from masses and radii. 

We first investigated the correlation of M+, R*, and logg 
with the effective temperature and with the stellar metallicity by 
using the Spearman correlation test. Results are given in Table 7. 


Article number, page 9 of 14 





























































A&A proofs: manuscript no. mdwarfs 




Fig. 12. Stellar mass (left panel), and radius (right panel), as a function of the effective temperature. Stars are plotted using different colours and 
symbols according to their metallicity. Several fits for fixed metallicity values are plotted: +0.15 (dashed line), +0.00 (solid line), -0.15 (dash- 
dotted line), and -0.30 (dotted line). The upper left panel shows the differences between the mass obtained from Equation 5 and those derived by 
using Henry & McCarthy (1993) relationship. The upper right panel shows the differences between the radius derived from Equation 6 and by 
using Equation 4. 


Table 7. Results from the Spearman’s correlation test and the MC 
method showing the dependence of the evolutionary parameters on the 
effective temperature, and on the stellar metallicity. 




Teff 



Spearman’s test 

MC method 


P 

prob. 

P 

Z-score 

Mass 

0.72 

2*10“'' 

0.47 

1.55 

Radius 

0.72 

2*10“® 

0.47 

1.55 

logg 

-0.73 

2*10“® 

-0.30 

-0.93 



[Fe/H] 



Spearman’s test 

MC method 


P 

prob. 

P 

Z-score 

Mass 

0.34 

0.02 

0.20 

0.60 

Radius 

0.34 

0.02 

0.20 

0.60 

logg 

-0.33 

0.02 

-0.12 

-0.40 


It can be seen that although the main dependence of the evo¬ 
lutionary parameters is on the effective temperature, they also 
show a moderate but significant dependence on the stellar metal¬ 
licity. 

We also evaluated the significance of the correlations by a 
bootstrapp Monte Carlo (MC) test plus a gaussian, random shift 
of each data-point within its error bars. For each pair of variables 
10000 random datasets were created, determining the coefficient 
of correlation, p, and its corresponding z-score each time. The 
tests were done using the code MCSpearman® by Curran (2014) 
which might be consulted for further details on this method. 
Mean values of p and z-score are given in Table 7. We note that 
the results from the MC method do not exclude a metallicity de¬ 
pendence but suggest that any possible correlation is relatively 
weak. 


® https://github.eom/PACurTan/MCSpearman/ 


A dependence of the radius on stellar metallicity is expected 
from model predictions (Baraffe et al. 1998; Dotter et al. 2008) 
however BO 12 find the interferometry-based radius rather insen¬ 
sitive to metallicity. Furthermore, while fitting M*, R*, and L* 
as a function of the effective temperature Mann et al. (2013b) 
find that adding the stellar metallicity as a parameter does not 
improve the fits. However, in a more recent work, Mann et al. 
(2015) do find a significant effect of the metallicity on the Tgff- 
R* relation. The authors point towards small sample sizes and 
a sparser sampling on [Fe/H] as the reasons why the effect of 
[Fe/H] was not noticed in their previous studies. We therefore 
performed two kinds of fit, one using only the effective temper¬ 
ature and another one adding the stellar metallicity as a parame¬ 
ter. The extra sum-of-squares F test (e.g. Lupton 1993) was used 
to test whether the addition of the metallicity to the functional 
form of the calibrations provides any improvement or not. The 
test returns a measure of the likelihood (p-value) that the simpler 
model (the one with fewer parameters) provides a better repre¬ 
sentation than the more complicated one. The resulting values'®, 
p-M* ~ 2x10“^, p-R* ~ 7x10“®, p-logg ~ 2x10“^, indicate 
that by including the metallicity there is improvement in the fits 
in line with Mann et al. (2015). The relationships we obtain are 
the following: 

M*(Mo) = -171.616 + 0.139 X T^s - 3.776 x lO^^T^j, 

+3.419 X lO^^r^jf + 0.382 x [FejH] (5) 
P*(Pq) = -159.857 + 0.130 x - 3.534 x lO^^T^j, 

+3.208 X 10“®r®3. + 0.347 x [FelH] (6) 
logg(cgi) = 174.462 - 0.138 x T^g + 3.728 x lO^^T^g: 

-3.376 X lO^^T^jf - 0.332 x [FejH] (7) 

The tests were performed using the MPFTEST IDL routine in¬ 
cluded in the MPFIT package (Markwardt 2009) and available at 
http://cow.physics.wisc.edu/ craigm/idl/idl.html 
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where the rms standard deviations of the calibrations are ctm* = 
0.02 Mq, ctr* = 0.02 Rq , and (Xiogg = 0.02 (cgs). The calibra¬ 
tions are valid for 3340 K < Teff < 3840 K, and -0.40 < [Fe/H] < 
-1-0.16 dex. Empirical relationships for stellar luminosity are not 
provided since they can be easily obtained from Teff and R* just 
applying Stefan-Boltzmann’s law. All these quantities (M*, R*, 
logg, and log(L*/Lo) ) for the stars analyzed in this work are 
provided in Table 8. 

Typical uncertainties are in the order of 13.1% for the stellar 
mass, 11.8% for the radius, 25% for luminosities, and 0.05 dex 
for logg. We note that these uncertainties are computed by tak¬ 
ing into account the cr of the corresponding calibration and the 
propagation of the errors in Teff and [Fe/H]. A word of caution 
should be given regarding the relative errors in mass since they 
tend to increase towards lower masses. Relative errors in mass 
might be larger than 20% for stars with M* < 0.35 Mq and reach 
up to more than 40% for the few stars with M* < 0.25 Mq. In a 
similar way, relative errors in radius can be larger than 20% for 
stars with R* < 0.35 R©. Relative errors in luminosities are also 
larger for low-luminosity stars, being significantly high (larger 
than 70%) for those stars with log(L*/Lo) < -2. We point as a 
possible explanation the fact that relative errors in masses ob¬ 
tained from Henry & McCarthy (1993) relationship tend to be 
larger at lower masses. 

The M* and R* versus temperature planes are shown in Fig¬ 
ure 12 where the stars are plotted with different colours accord¬ 
ing to their metallicities. It can be seen that for a given effec¬ 
tive temperature, larger stellar metallicities predict larger masses 
and radii. Regarding surface gravity, see Figure 13, the effect of 
metallicity tends to be the opposite with lower gravities in stars 
with higher metallicity content. 



Fig. 13. Logarithmic surface gravity as a function of the effective tem¬ 
perature. Stars are plotted using different colours and symbols according 
to their metallicity. Several fits for fixed metallicity values are plotted: 
+0.15 (dashed line), +0.00 (solid line), -0.15 (dash-dotted line), and - 
0.30 (dotted line). The upper panel shows the differences between the 
log g obtained with our calibrations (Equation 7) and those derived from 
masses and radius. Median errors in log g estimates are 0.18 dex (mass- 
radius derived values) and 0.05 dex (values from Equation 7). 


6. Summary 

The determination of accurate stellar parameters of low-mass 
stars is certainly a major topic in nowadays astrophysics. This 
is in part because of their advantages with respect to solar-type 
stars in the search for small, rocky, potentially habitable planets. 
This fact motivated us to develop a methodology to characterise 
early M dwarfs using the high-resolution spectra that are being 
obtained in the current radial velocity exoplanet programmes. 
We made use of ratios of features as a method to determine effec¬ 
tive temperatures, and combinations of features and temperature- 
sensitive ratios to determine metallicities. This technique largely 
applied to solar-type, subgiant, and giant stars had not been ex¬ 
tended before to the low-mass stars regime, probably because of 
the difficulty in identifying spectral lines in their optical spec¬ 
tra. We also provide empirical calibrations for masses, radii, and 
gravities as a function of effective temperature and metallicity. 
Our main results are as follows: 

- The behaviour of the EW of features was studied as a func¬ 
tion of the effective temperature and the metallicity. The re¬ 
sults show that for a significant fraction of the features, ~ 
50%, the EW shows a high anticorrelation with Teff, whilst 
the correlations between EW and metallicity are in general 
weaker. 

- Empirical calibrations for the effective temperature were ob¬ 
tained using stars with interferometric measurement of their 
radii from B012 as calibrators. 112 ratios of features sen¬ 
sitive to the temperature were calibrated providing effective 
temperatures with typical uncertainties of the order of 70 K. 

- In the same way 82 ratios of features were calibrated to de¬ 
rive spectral types. 

- Stellar metallicities were obtained from 696 combinations of 
EW of individual features and temperature-sensitive ratios, 
with estimated uncertainties in the range of 0.07-0.10 dex. 

- We made use of our technique to characterise 53 early M 
dwarfs which are currently being monitored in the HARPS 
exoplanet search programme. Photometric estimates of stel¬ 
lar mass, radius, and surface gravity were used to search for 
possible correlations of these parameters on Tgff and [Fe/H]. 

- We found stellar masses, radii, and surface gravities to have a 
moderate but statistically significant correlation with the stel¬ 
lar metallicity, in the sense that at a given effective temper¬ 
ature larger metallicities predict slightly larger masses and 
radii whereas, larger gravities are found in stars with lower 
metallicity content. 

Although high-resolution HARPS and HARPS-N optical 
spectra were used for this work, a similar methodology can be 
used to derive Teff and [Fe/H] for other instruments/spectral 
ranges. 
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Table 8 lists all the stars analyzed in this work. The table pro¬ 
vides: Star identifier (column 1); effective temperature in kelvin 
(column 2); spectral type (column 3); stellar metallicity in dex 
(column 4); logarithm of the surface gravity, logg, in cms“^ (col¬ 
umn 5); stellar mass in solar units (column 6); stellar radius in 
solar units (column 7); and stellar luminosity, log(L*/Lo) (col¬ 
umn 8). Each measured quantity is accompanied by its corre¬ 
sponding uncertainty. 
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Table 8. Stellar parameters for the stars studied in this work. 


Star 

(1) 

Teff 

(K) 

(2) 

Sp-Type 

(3) 

[Fe/H] 

(dex) 

(4) 

logg 

(cgs) 

(5) 

M* 

(Mo) 

(6) 

R* 

(Rg) 

(7) 

log(L*/Lo) 

(8) 

Gll 

3482 + 68 

M2.5 

-0.27 + 0.09 

4.91 + 0.06 

0.33 + 21.39% 

0.33 + 18.84% 

-1.834+ 38.49 % 

G187 

3562 + 68 

M2 

-0.16 + 0.09 

4.83 + 0.05 

0.43 + 13.22% 

0.42 + 11.99% 

-1.590 + 25.18 % 

Gll 76 

3603 + 68 

M2 

0.03 + 0.09 

4.75 + 0.04 

0.52 + 10.07 % 

0.51+9.31 % 

-1.409 + 20.08% 

G1191 

3587 + 68 

MO 

-0.39 + 0.09 

4.89 + 0.05 

0.35 + 15.26 % 

0.35 + 13.64% 

-1.728 + 28.31 % 

G1205 

3800 + 68 

M1.5 

0.00 + 0.09 

4.68 + 0.05 

0.60 + 9.35 % 

0.58 + 8.99 % 

-1.194+ 19.36% 

G1229 

3779 + 69 

Ml 

-0.10 + 0.09 

4.72 + 0.05 

0.55 + 9.81 % 

0.54 + 9.34 % 

-1.276 + 20.04% 

HIP31293 

3526 + 68 

M3 

-0.04 + 0.09 

4.81+0.05 

0.45 + 13.63 % 

0.44 + 12.39 % 

-1.564 + 25.95 % 

G1250B 

3557 + 68 

M2.5 

-0.13 + 0.09 

4.82 + 0.05 

0.44+ 13.08% 

0.43 + 11.88 % 

-1.576 + 24.97% 

G1273 

3342 + 69 

M4 

0.01 + 0.09 

4.97 + 0.11 

0.27 + 44.09 % 

0.28 + 37.91 % 

-2.055 + 76.27 % 

GJ2066 

3575 + 68 

M2 

-0.17 + 0.09 

4.82 + 0.05 

0.43 + 12.75 % 

0.43 + 11.59% 

-1.577 + 24.40% 

G1341 

3783 + 69 

M0.5 

-0.16 + 0.09 

4.74 + 0.05 

0.53 + 10.27 % 

0.52 + 9.75 % 

-1.305 + 20.81 % 

G1357 

3477 + 68 

M2.5 

-0.27 + 0.09 

4.92 + 0.06 

0.32 + 21.83 % 

0.33 + 19.21 % 

-1.844+ 39.21 % 

G1358 

3450 + 68 

M3 

0.04 + 0.09 

4.84 + 0.07 

0.42+ 18.71 % 

0.42 + 16.86 % 

-1.660+ 34.63 % 

G1367 

3559 + 68 

M2.5 

-0.06 + 0.09 

4.80 + 0.05 

0.46 + 12.22 % 

0.46+ 11.15% 

-1.525 + 23.59 % 

G1382 

3653 + 68 

M2 

-0.01 + 0.09 

4.74 + 0.04 

0.53 + 9.48 % 

0.52 + 8.82 % 

-1.372+ 19.15 % 

G1388 

3473 + 68 

M3.5 

0.12 + 0.10 

4.79 + 0.07 

0.47 + 15.76 % 

0.46 + 14.36 % 

-1.555 + 29.77% 

G1393 

3544 + 68 

M2 

-0.17 + 0.09 

4.84 + 0.05 

0.41 + 14.23 % 

0.41 + 12.86% 

-1.626 + 26.84% 

G1413.1 

3570 + 68 

M2 

-0.10 + 0.09 

4.80 + 0.05 

0.45 + 12.21 % 

0.45 + 11.14% 

-1.536 + 23.56 % 

G1433 

3618 + 68 

M1.5 

-0.13 + 0.09 

4.79 + 0.04 

0.47 + 11.02% 

0.46+ 10.13% 

-1.490 + 21.61 % 

G1438 

3647 + 68 

Ml 

-0.27 + 0.09 

4.83 + 0.04 

0.43 + 11.78% 

0.42 + 10.78 % 

-1.548 + 22.82% 

G1447 

3382 + 68 

M4 

-0.25 + 0.09 

5.00 + 0.09 

0.23 + 44.05 % 

0.24 + 36.99 % 

-2.156 + 74.41 % 

G1465 

3403 + 69 

M2.5 

-0.33 + 0.09 

5.01 + 0.09 

0.23 + 41.42% 

0.24 + 34.72 % 

-2.154 + 69.91 % 

G1479 

3476 + 68 

M3 

0.00 + 0.09 

4.83 + 0.06 

0.43 + 16.71 % 

0.42+ 15.10% 

-1.630+ 31.20% 

G1514 

3728 ± 68 

Ml 

-0.14 + 0.09 

4.76 + 0.04 

0.51+9.92% 

0.50 + 9.30 % 

-1.362 + 20.00% 

G1526 

3609 ± 68 

M2 

-0.10 + 0.09 

4.79 + 0.04 

0.47 + 10.99 % 

0.47 + 10.10% 

-1.483 + 21.57% 

G1536 

3685 ± 68 

Ml 

-0.08 + 0.09 

4.75 + 0.04 

0.52 + 9.66 % 

0.50 + 9.00 % 

-1.377+ 19.47% 

G1551 

3555 + 68 

M4.5 

-0.03 + 0.09 

4.79 + 0.05 

0.47 + 12.04 % 

0.46+ 11.01 % 

-1.511 + 23.29 % 

G1569A 

3608 + 68 

M2 

-0.02 + 0.09 

4.76 + 0.04 

0.50 + 10.32 % 

0.49 + 9.52 % 

-1.433 + 20.48 % 

G1581 

3419 + 68 

M3 

-0.20 + 0.09 

4.95 + 0.08 

0.29 + 29.73 % 

0.30 + 25.84 % 

-1.948 + 52.29% 

G1588 

3525 + 68 

M3 

0.00 + 0.09 

4.79 + 0.05 

0.46 + 13.21 % 

0.46 + 12.04 % 

-1.539 + 25.28 % 

G1618A 

3451 +68 

M3 

-0.10 + 0.09 

4.88 + 0.07 

0.37 + 21.30% 

0.37 + 18.96% 

-1.764+ 38.74% 

G1628 

3345 + 69 

M3.5 

-0.05 + 0.09 

4.99 + 0.11 

0.25 + 46.49 % 

0.26 + 39.59 % 

-2.104 + 79.61 % 

GJ644A 

3463 + 68 

M2.5 

0.08 + 0.10 

4.81 + 0.07 

0.45 + 17.19% 

0.44 + 15.59 % 

-1.603 + 32.17% 

G1667C* 

3572 + 68 

Ml 






G1674 

3484 + 68 

M2.5 

-0.20 + 0.09 

4.89 + 0.06 

0.36 + 19.49 % 

0.36+ 17.32% 

-1.767 + 35.52 % 

G1678.1A 

3815 + 69 

MO 

-0.09 + 0.09 

4.70 + 0.05 

0.58+ 10.13% 

0.56 + 9.73 % 

-1.222 + 20.77% 

G1680 

3585 + 68 

M2 

-0.12 + 0.09 

4.80 + 0.05 

0.45 + 11.84% 

0.45 + 10.82 % 

-1.527 + 22.93 % 

G1682 

3393 + 69 

M4 

0.02 + 0.09 

4.90 + 0.09 

0.35 + 28.14% 

0.35 + 24.92 % 

-1.834 + 50.50% 

G1686 

3677 + 68 

Ml 

-0.26 + 0.09 

4.81 + 0.04 

0.44+ 11.23% 

0.44 + 10.34 % 

-1.502 + 21.97% 

G1693 

3390 + 68 

M3 

-0.27 + 0.09 

5.00 + 0.09 

0.23 + 42.13% 

0.25 + 35.45 % 

-2.141 + 71.36% 

G1701 

3664 + 68 

Ml 

-0.19 + 0.09 

4.80 + 0.04 

0.46 + 10.75 % 

0.46 + 9.92 % 

-1.471 + 21.19% 

G1729 

3548 + 68 

M3.5 

-0.06 + 0.10 

4.80 + 0.05 

0.46+ 13.20% 

0.45 + 12.04 % 

-1.542 + 25.27% 

G1752A 

3551 +68 

M3 

0.02 + 0.09 

4.77 + 0.05 

0.49+ 11.78% 

0.48 + 10.79 % 

-1.485 + 22.91 % 

G1803* 

3628 + 68 

M0.5 






G1832 

3580 + 68 

M2 

-0.16 + 0.09 

4.82 + 0.05 

0.44 + 12.45 % 

0.43 + 11.33 % 

-1.562 + 23.91 % 

G1846 

3835 + 69 

MO 

-0.09 + 0.09 

4.69 + 0.05 

0.59 + 10.43 % 

0.57 + 10.07 % 

-1.194 + 21.39% 

G1849 

3486 + 69 

M3.5 

0.12 + 0.09 

4.78 + 0.06 

0.48 + 14.39 % 

0.47+ 13.13% 

-1.529 + 27.41 % 

G1876 

3357 + 68 

M4 

0.16 + 0.09 

4.90 + 0.10 

0.35 + 31.71 % 

0.35 + 28.14% 

-1.845 + 56.86 % 

G1877 

3428 + 68 

M3 

-0.11 + 0.09 

4.91 + 0.08 

0.34 + 25.11 % 

0.34 + 22.17% 

-1.838 + 45.06% 

G1880 

3736 + 68 

M1.5 

-0.01 + 0.09 

4.71 + 0.04 

0.57 + 9.04 % 

0.55 + 8.55 % 

-1.278 + 18.59 % 

G1887 

3712 + 68 

Ml 

-0.18 + 0.09 

4.78 + 0.04 

0.49 + 10.26 % 

0.48 + 9.57 % 

-1.406 + 20.50% 

G1908 

3570 + 68 

M1.5 

-0.27 + 0.09 

4.86 + 0.05 

0.39 + 14.23 % 

0.39 + 12.82 % 

-1.658 + 26.75 % 

LTT9759 

3581 +68 

M2.5 

0.07 + 0.09 

4.74 + 0.05 

0.52 + 10.33 % 

0.51+9.54% 

-1.413 + 20.54% 


Notes. * The star falls out of the range of applicability of our metallicity calibrations. 






